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A Landau model is used to study the phase behavior of the surface layer for magnetic and 
cholesteric liquid crystal systems that are at or near a Lifshitz point marking the boundary between 
modulated and homogeneous bulk phases. The model incorporates surface and bulk fields and 
includes a term in the free energy proportional to the square of the second derivative of the order 
parameter in addition to the usual term involving the square of the first derivative. In the limit of 
vanishing bulk field, three distinct types of surface ordering are possible: a wetting layer, a non-wet 
layer having a small deviation from bulk order, and a different non-wet layer with a large deviation 
from bulk order which decays non-monotonically as distance from the wall increases. In particular 
the large deviation non-wet layer is a feature of systems at the Lifshitz point and also those having 
only homogeneous bulk phases. 
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PACS numbers: 61.30.Cz; 64.60.Kw; 64.70.Md 
I. INTRODUCTION 

The interaction of a bulk system with a wall may give 
rise to a large variety surface phenomena, associated with 
the thermodynamic behavior of the surface layer adjacent 
to the wall. For example in ferromagnetic systems, when 
the interaction with the wall is such that it enhances lo- 
cal order it may happen that a surface transition takes 
place at temperatures above the critical temperature of 
the bulk. In such a transition the layers close to the wall 
become ordered although the bulk remains disordered. 
Depending on the nature of the interactions within the 
bulk and the interactions between the bulk and the wall, 
the system may exhibit phenomena such as wetting, crit- 
ical wetting, prewetting and other surface phase transi- 
tions. These phenomena have been extensively studied, 
both theoretically and experimentally in recent years (for 
a review see Q). 

A study of the global phase diagram for surface critical 
phenomena in ferromagnetic and other homogeneously 
ordered systems has been carried out by Nakanishi and 
Fisher || . In this study, a Landau phenomenological ap- 
proach has been applied and the phase diagram has been 
analyzed in the space of temperature, surface enhanced 
interactions, and bulk and surface ordering fields. For 
example it has been found that for finite positive surface 
field and no surface enhanced interactions, and in the 
limit of vanishingly small negative bulk field, the system 
exhibits a wetting transition as the temperature is varied 
below the bulk ordering temperature. At low temper- 
atures the surface field induces a local order in a layer 
of finite thickness I near the wall. However, at tempera- 
tures just below the bulk ordering temperature, the thick- 
ness surface layer is infinite, yielding a "wet" state. The 
two regimes are separated by a first order transition in 
which the thickness of the layer undergoes a discontinu- 



ous jump. This is known as the wetting transition. 

More recently surface phenomena in modulated systems 
have been considered. These systems are characterized 
by a periodic spatial variation of the order parameter 
in the bulk. Examples are magnetic spirals, cholesteric 
liquid crystals, amphiphilic systems, diblock copolymers 
and many others. In many cases the modulated phase 
is driven by a gradient-squared term with negative co- 
efficient in the Landau free energy. The system is then 
stabilized by terms quadratic in the second derivative. 
Studies of surface phenomena in such systems suggest 
that surface phase diagrams are rather rich, exhibiting 
novel surface states || and complicated surface struc- 
tures However the possible global phase diagrams 
of these systems have not been fully explored. 

Systems exhibiting a Lifshitz point may be considered 
as intermediate between ferromagnetic and modulated 
. In the Landau free energy of such systems, the 
coefficient of the gradient-squared term vanishes, making 
the quadratic term in the second derivatives the leading 
order interaction term. Surface phase diagrams of these 
systems have not been explored so far and it would be of 
interest to study them in some detail. 

In this paper we study the surface states and the sur- 
face phase diagram corresponding to a model of a Lifshitz 
point within the Landau approach. The phase diagram 
is studied in the space of temperature, bulk and surface 
fields. It is found that unlike the ferromagnetic case, 
these systems do not exhibit a wet phase in which the 
thickness of the surface layer diverges. It rather exhibits 
a transition from one surface state to another, as the 
temperature is varied, where both surface states have a 
finite thickness. 

We also consider the surface phase diagram of a fer- 
romagnetic system which is characterized by higher or- 
der interaction terms. Specifically we consider a Landau 
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free energy which includes terms quadratic in the gradi- 
ent and in the second derivatives of the order parame- 
ter, both with positive coefficient. As is well known, the 
quadratic term in the second derivatives does not affect 
the bulk phase diagram as long as the sign of its coef- 
ficient is positive. However we find, rather surprisingly, 
that although this term does not introduce any compet- 
ing with the gradient squared term, it affects the surface 
diagram in a profound way. In particular, we find that in 
addition to the usual wet and non-wet states which exist 
in the model of Nakanishi and Fisher, the model exhibits 
a second non-wet state with a distinct structure of the or- 
der parameter near the surface. Numerical studies yield 
the global phase diagram of the model. 

The paper is organized as follows: in Section II we 
review the model of Nakanishi and Fisher, and present 
analytic expressions for the location of the wetting and 
the critical prewetting points. Results of a numerical 
study of the surface phase diagram corresponding to the 
model with a Lifshitz point are presented in Section III. 
The surface states and the surface phase diagram of the 
generalized ferromagnetic model are discussed in Section 
IV. Finally a short summary is given in Section V. 

II. FERROMAGNETIC MODEL 

In this section we consider the surface phase diagram 
of a ferromagnetic system in the space of temperature, 
and bulk and surface ordering fields. The phase diagram 
exhibits wetting, prewetting and critical prewetting tran- 
sitions. The Landau phenomenological model of Nakan- 
ishi and Fisher is reviewed, and analytic expressions for 
the wetting and the critical prewetting points are given. 
In this model the ferromagnetic interaction is simply in- 
troduced by a gradient-squared term with positive co- 
efficient. Extensions of this model to other ferromag- 
netic systems with higher order ferromagnetic interac- 
tions, and to systems exhibiting Lifshitz points will be 
discussed in the following sections. 

Let 4>(x) be the scalar order parameter which corre- 
sponds to the ferromagnetic order. The wall is taken to 
be in the y — z plane, and the order parameter is assumed 
to depend only on the coordinate x perpendicular to the 
wall. The Landau free energy is given by 

F= [ dx{~W + \rcj) 2 + \^ + \(cj)'f) - h s c^ s (1) 
Jo 

where <j) s = </>(0) denotes the value of the order parameter 
at the wall and cf>' — dc/)/dx. The system is assumed to be 
of length L in the x direction. In calculating the surface 
free energy we will take the limit L — > oo. We have 
scaled the order parameter, the energy and the unit of 
length to simplify the coefficients, and so the bulk field h, 
the temperature r, and the surface field h s are rescaled 
variables. 
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FIG. 1. A schematic (r, h) phase diagram of the Nakan- 
ishi-Fisher model for h s > 0. The horizontal line represents 
the bulk coexistence line, which terminates at an ordinary 
critical point C. The curved line at negative bulk fields is the 
prewetting transition line on which the two types of surface 
solutions coexist. This line intersects the h = axis at the 
wetting point WT which separates the dry phase existing at 
large negative r, known as the prewet phase, PW from the wet 
phase, W, which exists close to the critical point. The prewet- 
ting transition line terminates at the critical prewetting point 
CP, where the two surface solutions become identical. The 
width I of the surface layer of the wet and the prewet phases 
is finite for non-vanishing bulk field. However / diverges in 
the limit h — > CP in the wet phase. 

The (r, h) phase diagram of this model for non- 
vanishing surface field h s > is given schematically in 
Figure 1 j^]. The order parameter away from the wall 
is not affected by the surface field. For a negative bulk 
field, h < 0, it approaches a negative value character- 
istic of the bulk. On the other hand the surface field 
enhances the order within a layer of thickness I. The 
thickness of this layer undergoes a discontinuous change 
along the prewetting transition line in the Figure. For 
h = this becomes the first order wetting transition, 
WT. The line ends at some critical point CP known as 
the critical prewetting point. The low temperature sur- 
face state, existing to the left of the line is the prewet 
state, PW. It is characterized by a surface layer with 
a finite width. The state to the right of the line is the 
wet state, W, in which the width of the surface layer 
diverges in the limit of vanishing bulk field. In the fol- 
lowing we analyze the Landau free energy (|l|) and obtain 
analytic expressions for the wetting WT and the criti- 
cal prewetting CP points. The Euler-Lagrange equation 
corresponding to the free energy [l] is 

<t>" + h - rep - 3 = (2) 

with the boundary condition at x = 

hs = ~K (3) 
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FIG. 2. Schematic order parameter profiles for h < and 
h s > 0, (a) at large negative r and (b) close to the bulk critical 
line. In the limit h —* 0~ , the width of the surface layer of 
solution (b) diverges, leading to a wet state. 



We are interested in calculating the order parameter pro- 
file and free energy for negative bulk field h < and pos- 
itive surface field h s > 0. We thus expect that for large 
x, the order parameter approaches the bulk value —<pB 
where 4>b > satisfies 



h + rcj>B + 







(4) 



Multiplying Eq. 
grated to yield 



by <[)' this equation may be inte- 
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where C is a constant. This constant may be evaluated 
by noting that at large x the order parameter asymptot- 
ically approaches —<f>B- Thus 



c = W B + u% + h<t> 



(6) 



Using this result, the first integral of the Euler equation 
(0) becomes 
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where the (— ) sign in the right hand side is taken since for 
the choice of the bulk and surface fields in this problem 
the order parameter is expected to decrease with x. 

In order to locate the wetting point WT we take the 
limit h -> CT in Eq. (0) 
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In order to evaluate the surface free energy, F s , associated 
with the local surface order we note that the free energy 
density of the bulk state is given by C. Thus F s = F — 
CL. Using (0), (g) and (g) one obtains 



or 



F s = F-CL 



F* = 



dx((p') 
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(9) 
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We proceed by evaluating the order parameter profile 
obtained from Eq. (g). Two distinct types of profiles 
are found: 4>i(x) for negative and large r and 4>2(x) for 
negative and small r, close to the bulk critical point. 
These profiles are schematically given in Figure 2. For 
large r, the surface field does not affect the local surface 
order in a substantial way and thus the surface order 
parameter <p s i remains close to the bulk value, satisfying 
—(pB < 4>si < 4>b- Integrating (||) one finds the surface 
free energy for this type of solution 
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where the surface order parameter, cf> s i, is determined by 
the boundary equation 
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On the other hand for small r the local order is highly 
susceptible to the local ordering field and one obtains 
an order parameter profile 4>2{x) which at the surface 
is considerably different from the bulk value, satisfying 
0s2 > 4>b- Integrating (||) for this solution one obtains 
the surface free energy 
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with the surface order parameter satisfying 



(13) 
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At the wetting transition the two types of solutions have 
the same free energy. To find the transition point we 
define y\ = 4> s i/4>Bi 2/2 = 4 > s2/4'b- At the coexistence 
point one has 
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In addition, the boundary condition equations (12) and 
( |l4| ) impose the following relation between y\ and yi : 

vl +2/2=2 (16) 

To solve Eqs. © and (|l|) we use the substitutions 
y\ = 1 — u and y 2 = 1 + u where u — h s ^/2/ <f) 2 E 
( |l5| ) is then readily solved yielding u 2 = \/l2 
wetting transition thus takes place at 



B . Equation 
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To locate the critical prewetting CP point we consider 
the boundary equation (|^). Combining it with the ex- 
pression for the order parameter derivative (Q) it may be 
written as 

G(cj> s ) = <£ + 2r<t>l - 4hcj> s + (<f>% - 2hcp B - 2h 2 s ) = 

(18) 

This equation determines the surface order parameter <f> s 
for given r,h and h s . At the critical prewetting point 
the two solutions of this equation, which correspond 
to the two coexisting states on the prewetting line be- 
come identical. The conditions for this to take place are 
dGJ d<p s — d 2 G/d<fi 2 = 0. These equations together with 
( plf ) yield the critical prewetting point 
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h = -2(^)3/4^3/2 
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It is of interest to explore the general validity of this 
global phase diagram by considering Landau free energies 
with different non-linear terms. To this end we stud- 
ied the phase diagram of a Landau model with piecewise 
parabolic potential 
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dx(-h(f> + f(<f>) 
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where 
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and the parameters a and are dependent on r. The 
analysis presented above for the Nakanishi-Fisher model 
may be extended to study the phase diagram of the Lan- 
dau free energy ( |2C| ) . It is found that the (r, h) phase 
diagram of the two models exhibit the same qualitative 
features. The wetting transition takes place at 
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(22) 



while the critical prewetting point is found to be located 
at 
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(23) 



III. LIFSHITZ POINT MODEL 

In this section we study the surface phase diagram cor- 
responding to a model of a Lifshitz point in the presence 
of a surface ordering field. We consider the Landau free 
energy 

F = -h s 4> s 

+ [ dx{~h(j) + \r<p 2 + \^ + \v{fif + i(0") 2 ) (24) 
Jo 

with v = 0. This model for v > will be considered 
in the next section. The Euler-Lagrange equation corre- 
sponding to this free energy is 



- v<j>" - h + r(f> + (f) 3 =0 



with the boundary equations at 
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We have integrated numerically Eq. (|25j) with the 
boundary conditions (26) for v = 0, corresponding to the 
case of a Lifshitz point. We find that unlike the ferromag- 
netic case, the model does not exhibit a wet phase. It 
rather exhibits two distinct phases (a prewet state, PW, 
and a state, S, with surface enhanced order) both char- 
acterized by surface states with a finite width. Here, in 
analogy with the phase diagram of the Nakanishi-Fisher 
model we keep referring to the low temperature state as 
the prewet phase, although the phase diagram does not 
exhibit a wet state at all. The resulting (r, h) phase di- 
agram for h s = 0.1 is given in Fig. 3. The diagram 
displays a first order transition line separating the two 
surface phases PW and S. The line terminates at a sur- 
face critical point SC. This line intersects the h — axis 
at the point PS. 

Representative order parameter profiles in the two non- 
wet phases are given in Fig. 4 for small bulk field h. At 
low temperatures, r <C — 1, the surface field introduces 
only a weak local order near the surface, and the order 
parameter decays monotonically to its bulk value as one 
moves away from the wall. This is very similar to the 
low temperature phase of the ferromagnetic case. How- 
ever at higher temperatures, just below the bulk critical 
point, the order parameter becomes highly susceptible 
to the local surface field, the surface order parameter is 
much larger than the bulk value, and it decays in a non- 
monotonic way to the bulk value away from the wall. The 
width of the surface layer remains finite even in the limit 
h — > 0~. This is in a sharp contrast with the ferromag- 
netic case where the width diverges in this limit, leading 
to a wet state. 



IV. EXTENDED FERROMAGNETIC MODEL 

In this section we consider the extended ferromagnetic 
model (p3) with v = 1. We restrict this study to negative 
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FIG. 3. The (r, h) phase diagram corresponding to the Lif- 
shitz point model for h s = 0.1. The phase diagram displays 
two non-wet states: a prewet state, PW, and a state with 
surface enhanced order, S. The two phases are separated by 
a first order line. This line intersects the h — axis at the 
point PS. It terminates for finite ftata surface critical point 
SC. 
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FIG. 4. Characteristic profiles of the order parameter 
of the two non-wet states PW and S. The profiles 
are given at a point on the coexistence line of the two 
states, close to the point PS of Fig. 3 where h — > CP, 
r = -0.52008, h s = 0.1, h = -10~ 6 . 
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FIG. 5. The (r,h s ) phase diagram of the extended ferro- 
magnetic model in the limit h — > - . It displays three distinct 
states: a wet state, W, and two non-wet states PW and S. 
The state S is characterized by surface enhanced order (see 
Fig. 6). The phases are separated by three first order line 
which intersect at a triple point TP. The bulk critical point 
is denoted by C. 

bulk fields in the limit h — > CP and evaluate the (r,h s ) 
phase diagram. 

For small surface fields the model is found to yield sim- 
ilar surface phenomena as the model of Nakanishi and 
Fisher. It exhibits a wet phase, W, at temperatures be- 
low the bulk critical point and a prewet state, PW at 
low temperatures. The two states are separated by a 
first order wetting transition. 

However the phase diagram becomes rather different 
for large h s . Here, in addition to the wet state existing at 
high temperatures, one finds two distinct surface phases 
at lower temperatures: a phase with a surface enhanced 
order, S, and a prewet phase, PW. The two phases are 
separated from each other by a first order transition. The 
resulting (r, h s ) phase diagram is given in Fig. 5. We 
also display some characteristic order parameter profiles 
of the three phases. Fig. 6 gives the profiles at a point 
on the PW — S coexistence line, while the profiles at a 
point on the PW — W line are given in Fig. 7. 

V. SUMMARY 

Using a general Landau model, we have studied the 
surface phase diagrams that result from the effect of a 
wall bounding a semi-infinite sample which exhibits ho- 
mogeneous bulk phases. The model is applicable to a 
wide variety of systems including magnetic materials and 
cholesteric liquid crystals. 

Choosing coefficients to simplify the model to the ferro- 
magnetic model studied by Nakanishi and Fisher, we ob- 
tained analytic expressions for the temperature and bulk 
field co-ordinates of the wetting transition point and the 
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FIG. 6. Characteristic profiles of the order parameter at a 
point along the PW — S coexistence line for h — > CP. The 
profiles are given for r = -0.2, h s = 0.09967 and h = -10~ 10 . 
It is evident that while in the S state the surface order is 
highly susceptible to the surface field, the surface order in the 
PW state is affected only mildly by this field. 

critical prewetting point as functions of the surface field. 
We also demonstrated that the general surface phase di- 
agram is not highly sensitive to the precise form of the 
non-linear terms. 

In contrast, we found that a system at the Lifshitz 
point is similar to the case of modulated bulk phases 
where a wetting layer does not form. Instead, a non-wet 
surface state forms which decays non-monotonically. The 
wetting transition is replaced by a transition between the 
two types of non-wet states and the critical prewetting 
point becomes a surface critical point. 

In the extended ferromagnetic model we examined only 
the case of vanishing bulk field and found that all three 
types of surface layers develop if the surface field is not 
too weak. As temperature is reduced from the bulk crit- 
ical point, a wetting transition occurs, followed at even 
lower temperatures by a transition from a highly devi- 
ated to a weakly deviated non-wet layer. Unfortunately, 
we are not aware of any experimental observations of such 
a transition so far. This is, however, not surprising be- 
cause it is a subtle change that has not been expressly 
looked for. 
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FIG. 7. Characteristic profiles of the order parameter at 
a point along the PW — W coexistence line for h — > 0~ . 
The profiles are given for r — —0.05, h a = 0.02433 and 
h — — 10~ 10 . Note that in this figure the flat part of the 
W state has been truncated, as the width of the surface layer 
of this state diverges in the limit of vanishing bulk field. 
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